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Abstract 

The covariant Maxwell equations are derived from the continuity equation for the electric charge. 
This result provides an axiomatic approach to Maxwell's equations in which charge conservation 
is emphasized as the fundamental axiom underlying these equations. 
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1. Introduction 



In classical electrodynamics the conservation of electric charge, which is expressed through 
the continuity equation 

v.j + £-o, w 

is not usually considered to be an independent assumption but a formal consequence of 
Maxwell's equations. The familiar argument is that the divergence of the Ampere-Maxwell 
law together with the Gauss law yield equation (1). Regarding equation (1) Griffiths has 
written [1]: "This is, of course, the continuity equation — the precise mathematical statement 
of local conservation of charge. As I indicated earlier, it can be derived from Maxwell's 
equations. Conservation of charge is not an independent assumption, but a consequence of 
the laws of electrodynamics." The argument is certainly more elegant in spacetime where 
the covariant form of the continuity equation reads d u J u = 0. By applying the operator d v 
to the inhomogeneous Maxwell equations (in Gaussian units) d^F^ = (4ir/c)J u we obtain 
d u J u = because of the identity d u d^F^ u = 0. 

Nevertheless, we have recently questioned [2] the traditional idea that charge conservation 
is a consequence of Maxwell's equations. In most textbooks the Ampere-Maxwell law, 
one of the Maxwell equations, is introduced by assuming the validity of the continuity 
equation. This fact involves a circular reasoning: to obtain the Ampere-Maxwell law we 
assume charge conservation and therefore this conservation law should not be considered 
as a consequence of Maxwell's equations. Furthermore, from a historical point of view 
charge conservation was the ingredient which allowed Maxwell to generalize the quasi-static 
Ampere's law: V x B = (4n/c)J to the Ampere-Maxwell law: V x B - (l/c)dE/dt = 
(47r/c)J. Accordingly, charge conservation should be interpreted as an axiom of Maxwell's 
equations rather than a consequence of them. 

We have also demonstrated [2] an existence theorem for two retarded fields, which says 
that given localized time-dependent scalar and vector sources satisfying the continuity equa- 
tion there exist two retarded vector fields that satisfy four field equations. When the sources 
are identified with the usual electromagnetic charge and current densities and appropriate 
choice of constants made, the retarded vector fields become the generalized Coulomb and 
Biot-Savart laws [3,4] and their associated field equations become Maxwell's equations. We 
have concluded that charge conservation is the fundamental axiom underlying Maxwell's 



equations. Our derivation of Maxwell's equations has been commented on by Jefimenko 
[5,6] and Kapuscik [7,8]. Colussi and Wickramasekara [9] have stressed the result that 
Maxwell's equations, which are Lorentz-invariant, have been obtained from the continuity 
equation, which is Galilei-invariant according to these authors. 

An advantage of the derivation of Maxwell's equations presented in reference 2 is that it 
naturally introduces the time-dependent extensions of the Coulomb and Biot-Savart laws in 
the form given by Jefimenko [3,4]. However, a practical disadvantage of such a derivation is 
that it involves unfamiliar and large identities involving derivatives of retarded quantities. 

In this paper we extend our three-dimensional derivation of Maxwell's equations [2] to the 
four-dimensional spacetime. We show how the covariant form of Maxwell's equations can be 
obtained from the continuity equation for the electric charge. The existence theorem for two 
retarded fields demonstrated in reference 2 is now formulated for a retarded antisymmetric 
field tensor in spacetime. The existence theorem in spacetime is then applied to the usual 
electromagnetic four-current and thus the electromagnetic field together with the covariant 
Maxwell equations are obtained. This four-dimensional derivation of Maxwell's equations 
is considerably shorter than that formulated in the three-dimensional space. Readers well 
versed in the mathematics of the spacetime of special relativity could find more simple 
and elegant the four- dimensional derivation because it does not involve large operations 
on explicit retarded quantities but instead it involves simplified tensor operations in which 
retardation is implicit in the retarded Green function of the free-space wave equation. 

2. The covariant Maxwell equations 

Greek indices //, u, k . . . run from to 3; Latin indices k, . . . run from 1 to 3; x — 
x 11 = (ct, x) is the field point and x' = x ,fl = (ct' : x') the source point; the signature of 
the metric of spacetime is (+, — , — , — ); £^ a P i s the totally antisymmetric four- dimensional 
tensor with e 0123 = 1 and e' l ^ k is the totally antisymmetric three-dimensional tensor with 
£ i23 _ i Summation convention on repeated indices is adopted. 

The covariant form of Maxwell's equations can be written in Gaussian units as: 

4:71 

Q^V = _JU Q* W V = Q ? 

where F» v is the electromagnetic field tensor and *F^ U = (l/2)e^ Ka F Ka is the dual of F^ v . 
The antisymmetric tensor F^ v is defined by its components F l ° = (E) 1 and i™ = — e^ k {Q)k, 
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and therefore the components of *F^ V are *F t0 = (B)* and *F lj = s ljk (E) k . The source of the 
tensor F^ v is the electric four-current: J u = (cp, J), which is a conserved quantity 

d v r = o. (3) 

As above mentioned, equation (3) seems to be a formal consequence of equations (2) because 
of the identity dyd^F^ = 0, that is, equations (2) imply equation (3). We will investigate 
now the converse implication: To what extent equation (3) implies equations (2)? 

3. An existence theorem in spacetime 

Let us formulate the following existence theorem: Given the localized four-vector J v 
satisfying the continuity equation 

d v J v = 0, (4) 
there exists the antisymmetric tensor field 

= J d 4 x'G(d'»J u - d' u J»), (5) 

that satisfies the field equations 

Q^ = J u and 9^1^ = 0, (6) 

where G(x, x') = S(t' — t + R/c)/(AnR) with R = |x — x'| is the retarded Green function of 
the free-space wave equation and *J r ^ v = (l/2)e MiyQ/3 jF Q , ( g is the dual of J r ^ v ' . The function 
G satisfies the free-space wave equation & &^G = 5(x — x'), where 5(x — x') is the four- 
dimensional delta function. The integral in equation (5) is over all spacetime. We note that 
G is not an explicit Lorentz invariant object. This objection can be avoided by replacing G 
by the less-known retarded invariant form D(x, x') = 9(x — x' )S[(x — x') 2 }/ (2tt), where 9 is 
the theta function [10]. 

The proof of the theorem follows from the tensor identity 

d^G{d'»J u - d' u Jv)-d'^Gd w J^ - J v d^G)=J v 5(x - x'). (7) 

This identity is of general validity for a four-vector J v that satisfies equation (4) but oth- 
erwise is arbitrary. In particular, equation (7) is valid when J v is localized into a finite 
region of spacetime, property that guaranties that surface integrals in spacetime involving 
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J v vanish at infinity. Equation (7) can be obtained from Eq. (4) as follows. Assume that 
J v satisfies the continuity equation evaluated at the source point: 

V = 0. (8) 

We apply the operator —Gd' v to equation (8) and add the term ^d'^d'^G to both sides: 

-Gd^d'^ + T&^G = T&^G. (9) 

The left-hand side of equation (9) can be written as 

-^G{&^J V - d w J») - ^{Gd' u J» - J v d'»G) = J v d'^G. (10) 

When —d'^G = d^G is used in the first term of the left-hand side of equation (10) and 
d'^d'^G = 5(x — x') in its right-hand side, we finally obtain equation (7). 
We integrate equation (7) over all spacetime 

where <9 M has been extracted from the first integral of the left-hand side. The second integral 
on this side can be transformed into a surface integral which vanishes at infinity because of 
J v is assumed to be a localized quantity. The term on the right-hand side can be integrated 
yielding J v . Thus equation (11) reduces to 

d„J d 4 x'G{d , »J u - d' v J») = J v . (12) 

This equation shows the existence of the antisymmetric tensor 

= J d 4 x'G{d , »J v -&"J% (13) 

in terms of which equation (12) takes the compact (inhomogeneous) form: 

V 4 " = r. (14) 

The dual of T^ v can be written as 

* T »v = e »vK\ J d * x > Gd > K j x . ( 15 ) 

We apply the operator to equation (15) and perform an integration by parts to obtain 
d;^ v = e» UKX J d 4 x'Gd^d' K J x -J d 4 x'd'^ UKX Gd' K J x ). (16) 
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The first integral vanishes because of the symmetry of the indices ji and k in d'^& K and 
the antisymmetry of such indices in s^ ukX . The second integral can be transformed into 
a surface integral which vanishes at infinity. Then equation (16) reduces to the compact 
(homogeneous) form 

d*?^ = 0. (17) 

Once equations (13), (14), and (17) have been obtained, the theorem has been demonstrated. 
Notice that this theorem applies to any conserved four-vector which does not necessarily 
belong to the electromagnetic theory. 

4. Deriving Maxwell's equations 

Let us apply the existence theorem to the four- vector = (3J^, where J M = (cp, J) is 
the usual electric four-current with p and J being the charge and current densities and (5 is 
a constant. Then there exists the tensor field 

S r = pfMG i *-r-,rj>), (is, 

which satisfies the field equations 

d,g^ = pr and o;g^ = o, (19) 

where *Q^ V can be written as 

= f3e» UKX J d 4 x'Gd' K J x . (20) 

The constant (3 belongs to the a/37-system which allows us to write electromagnetic expres- 
sions in a way that is independent of units [11]. The constants a, (3 and 7 satisfy the relation 
a = P^c 2 with c being the speed of light in vacuum [11]. If we let a = l/e ,P = /i , and 
7 = 1, we obtain expressions in SI units, and if we let a = 47r, (3 = 4ir/c, and 7 = 1/c, we 
obtain expressions in Gaussian units. 

To verify that equations (19) are already the covariant Maxwell equations, that is, that 
the tensor identifies with the well-known electromagnetic field tensor F^ u , we proceed as 
follows. Let us label the components of the tensor Q^ v in equation (18) as Q M = (/3c/a:)(X) 1 
and Q l * = —e z: > k (Y) k , where X and Y are two time-dependent vector fields. The components 
of the dual tensor can be obtained from those of by making the dual changes 
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((3c/a)(X.y -> (Y)< and (Y) k -> -[l/( 7 c)](X) fc . Notice that /3c/a = l/( 7 c). It follows that 
*gM = (Y) j and = [l/( 7 c)]5^(X) fe . Our plan is now to identify the vectors X and Y 
with the electric and magnetic fields E and B. With this purpose in mind, let us consider 
the iO component of equation (18) which can be written as 



a 

47T 



d 4 x'G(d H j u -d'"j l ), 



-/ ld ^ d t'G\-iy' P r-- 2 [- 



c 2 V di- 



et 



4tt / R 



1 dJ 



(21) 



ret 



where the square brackets [ ] ret indicate that the enclosed quantity is to be evaluated at the 
retarded time t' — t — R/c. The right-hand side of Eq. (21) is the same as the right-hand 
side of the i-component of the retarded electric field 



(E)' = ±[ dV± 
K J Air J R 



p c 2 dt' 



(22) 



ret 



Comparison between equations (21) and (22) yields the identification X = E. We consider 
now the iO component of equation (20) which can be expressed as 



e l0a/3 J d A x'Gd' a Jp, 



47T 



(23) 



The right-hand side of equation (23) is the same that the right-hand side of the i-component 
of the retarded magnetic field 



(B) 



ref 



(24) 



Comparison between equations (22) and (23) yields the identification Y = B. Then the 
tensor Q^ v naturally identifies with the electromagnetic field tensor 



(25) 



which satisfies the Maxwell equations 



d^Fi" = /3J" and d*F^ = 0, 



(26) 
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where *F^ V can be written as 



= pe^ KX J d 4 x'Gd' K J x - (27) 

The components of F^ v are given by 

F i0 = —(E)* and F ij = -e i]k (B) k , (28) 
a 

where (E) 1 and (B) fc are the components of the electric and magnetic fields. It follows that 

*F i0 = (BY and *F ij = —e ijk (E) k . (29) 

7c 

Therefore, the covariant form of Maxwell's equations have been obtained as an application 
of the existence theorem formulated here. 

For completeness, we can verify that Eqs. (26) yield the three-dimensional form of 
Maxwell's equations in a/37-units. By using equations (28), (29) and <9 M = [(l/c)d/dt, V], 
we can write the four- vectors [11]: 

fl^=(*:V.E,VxB-2f\ (30) 
\a a at J 

a ^=(v.B,-lvxE-If). (31) 

From equations (26), (30), (31) and J u = (cp, J), we obtain the three-dimensional form of 
Maxwell's equations in af3^ units for sources in vacuum [11]: 

V • E = ap, (32) 

V • B = 0, (33) 

VxE + 1 — = 0, (34) 

VxB-^ = |JJ. (35) 
a at 

In particular, if we let a — 1/eo, f3 — Po, and 7 = 1, then we obtain Maxwell's equations in 
SI units, 

5. Discussion 



In section 4 we have established the identification Q^ v = F^ u on the basis that equations 
(22) and (24) are already known. But suppose that we do not know these equations but 
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only the experimental time-independent Coulomb and Biot-Savart laws. In this case the 
time-independent limit of the iO component of equation (18) can be written as 

W = £ / (36) 

where R = R/i2. From the experimental law / = g E for the force / acting on a charge q 
we know that there exists the electrostatic field E with components 

(E)'^/ d v(§> (37) 

Comparison between equations (36) and (37) yields the identification X = E in the time- 
independent regime. We consider now the time-independent limit of the iO component of 
equation (20) which can be expressed as 

W = £/«<•<'*!)' (38) 

From the experimental law / = 7 J d 3 x'J x B for the force / acting on a current J we 
know that there exists the magnetostatic field B given by 

( B >' = £/ dV ( JX :§)' < 39 > 

Comparison between equations (38) and (39) leads to the identification Y = B in the 
time-independent regime. 

Therefore, by requiring that the time-independent limit of equations (18) and (20) be 
consistent with the experimental Coulomb and Biot-Savart laws, we obtain the identifications 
X = E and Y = B in the time-independent regime. Such identifications naturally extend 
to the time-dependent regime of the theory and so we conclude again the identification 
gnu _ pnu _ us imagine for a moment that there is a physicist that only knows the 
experimental Coulomb and Biot-Savart laws. He then can infer the electromagnetic field 
together with Maxwell's equations by assuming the validity of the continuity equation. The 
physicist would correctly conclude that charge conservation is the fundamental postulate 
underlying Maxwell's equations. 

It can be argued that the derivation of the covariant form of Maxwell's equations from the 
continuity equation involves other implicit physical postulates like retardation (causality) 
which is implicit in the use of the retarded Green function G. If instead of the Minkowski 
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spacetime with metric signature (+, — , — , — ) together with its associated retarded Green 
function G = 5{t' — t + R/c}/(4irR), we consider, for example, an Euclidean four-space 
[12-14] with metric signature (+, +, +, +) together with its associated retarded-imaginary 
Green function [14] Gj = 5{t' — t + R/(ic)}/(4irR), then we obtain a set of Maxwell-like 
field equations in the Euclidean four-space [12-14]. This means that the continuity equation 
can also imply other electromagnetic theories different from that of Maxwell. This is an 
expected result because charge conservation is independent of the signature of the metric of 
spacetime. The independence of the signature can be verified in the equation 

d a j a = d a j a = v • J + ^ = dj a = d a j a , (40) 

where d a = [(l/c)d/dt,V\, d a = [(l/c)d/dt, -V], J a = (cp, - J) and J a = (cp, J) are 
objects in the Minkowski spacetime and d a = [(l/c)d/dt,'V], d a = [(l/c)d/dt, V], J a = 
(cp, J) and J a = (cp, J) are objects in the Euclidean spacetime (there is no difference be- 
tween covariance and contravariance in this spacetime). Hehl and Obukhov [15] have pointed 
out that charge conservation is metric-independent because it is based on a counting proce- 
dure for elementary charges. Therefore if one first postulates the validity of the continuity 
equation in any four-space, the particular use of the Minkowski spacetime (together with its 
associated retarded Green function G) would not really be a new postulate, but only just a 
particular application of the initial postulate. 

6. The subtle relation between sources and fields: Is this an egg-hen problem? 

A referee has pointed out that dp/dt and J can only be produced by electric and magnetic 
fields (take Ohm's law for instance J = crE). Therefore a circular process seems to be 
unavoidable in electromagnetism: p and J imply E and B which in turn imply new p\ and 
3 1, and so on. Because of this circular characteristic, it is not clear if E and B (satisfying 
Maxwell's equations) are a consequence of p and J (satisfying the continuity equation) or vice 
versa. According to the referee it seems a matter of taste to say which one is a consequence 
of the other. In other words: From referee's comment we could conclude that the connection 
between sources and fields is a little bit like the egg-hen problem: Who was first? 

Let us to briefly discuss the very deep relation between sources and fields in the context 
of the electromagnetic theory. It is now a common place to say that charges and currents are 
the causes of the electric and magnetic fields, but this was not so in the 19th century [16,17]. 
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Two rival point of views on what is now known as electromagnetic theory coexisted through 
most of the 19th century. On one hand, the field point of view followed by Faraday, Maxwell, 
Thomson, Fitzgerald, and Larmor who considered the fields (and their "observable" lines of 
force) to be the primary objects of the theory. They viewed the charges as manifestations 
of the terminal points of lines of force having no independent or substantial existence [17]. 
On the other hand, the charge-interaction point of view followed by Ampere, Neumann and 
Weber who considered the electromagnetic phenomena as originating in the interaction of 
charges which was regarded by some as mediated by an ethereal continuum and by others as 
a direct action at a distance. In the charge-interaction view, charges and currents were the 
primary objects of the theory and the causes of the forces. Modern field theory of the 20th 
century resulted from a combination of these two point of views. Charges were regarded 
as sources of the electromagnetic forces (this was the legacy of the charge-interaction point 
of view) which were considered as propagating through the Maxwellian field (this was the 
legacy of the field point of view). 

The existence of a causal relation between sources and fields is now universally accepted 
in electromagnetism: p and J (causes) produce the fields E and B (effects). This relation 
is clearly expressed by the retarded solutions of Maxwell's equations in the form of time- 
dependent extensions of the Coulomb and Biot-Savart laws (in Gaussian units) [3]: 



According to these equations the values of E and B at the instant t are determined by the 
values of p and J at the retarded (previous) time t — R/c. This result clearly agrees with 
the principle of causality which roughly speaking says that the causes (in this case p and J) 
precede in time to their effects (in this case E and B). Traditionally we place the causes on 
the right-hand side of equations and the effects on their left-hand side. 

The principle of causality provides then a reasonable argument to elucidate what are the 
causes and what are the effects in electromagnetism. Who were first? By looking at Eqs. (41) 
and (42), our answer is straightforward: p and J occurred first and therefore they should be 
considered the causes of E and B. This interpretation would seem incontrovertible — and 
so the supposed egg-hen problem involved in the source-field relation would be an illusion- 




's dtj 
J(x',f-fl/c)xR d_ 



p(x',t-R/c)R d 
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except because there exist also the so-called advanced solutions of Maxwell's equations, 
which can be written as 



/ 



A ,p(x-, f y/ c) R_| | rf ,^, t +fl/ c) R + j(^/ c ; ) = E(xi) (43) 

., ,J(x.',t+R/c)xR d f 3 ,J(x.',t+R/c) x R 



dV 1 g ; — a7 dV • ^ =B < x -» < 44) 

According to these equations the values of E and B at the instant t are related with the 
values of p and J at the advanced (posterior) time t+R/c. We again invoke the principle of 
causality to determine what are the causes and the effects in electromagnetism. Who were 
first? By looking at Eqs. (43) and (44) we note that E and B occurred first and therefore 
they should be considered the causes of p and J. Furthermore, if we introduce the variable 
r = t + R/c then Eqs. (43) and (44) can alternatively be written as 



/ 



/ 



3 ,J(x',t)xR d f 3 ,J(x',r)xR 



These equations say that the values of p and J at the instant r are related with the values 
of E and B at the previous time r—R/c. In the light of the causality principle, the advanced 
solutions of Maxwell's equations are consistent with the idea that E and B are the causes 
of p and J, but the retarded solutions of these same equations are consistent with the idea 
that p and J are the causes of E and B. As above mentioned, the subtle relation between 
sources and fields in electromagnetism is a little bit like the egg-hen problem. 

It can be argued, however, that the above analysis is purely formal and that other physical 
considerations allow us to identify the causes and the effects in electromagnetism. Suppose 
we have an electron at rest which is associated with an electrostatic field. In this case, 
we have neither retarded nor advanced times and so we cannot clearly distinguish between 
causes and effects. Suppose now that the electron is accelerated by mechanical forces. The 
resulting electric current obeys the continuity equation and the reason the charge moves 
has nothing to do with electromagnetism. Experience shows that the electron radiates, 
according to Maxwell's equations, and that this radiation carries of electromagnetic energy, 
which must come at the expense of the particle's mechanical energy. The radiated energy 
detaches from the electron, propagates off to infinity never coming back and we can physically 
detect this electromagnetic energy. In this dynamical case the radiation fields clearly arise 
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as a consequence of the motion of charges which was originated by mechanical and not 
electromagnetic means. This simple example shows that the continuity equation is valid 
outside the domain of Maxwell's equations and therefore it should be regarded as a stronger 
axiom. 

7. Concluding remarks 

R. G. Brown wrote [18]: "... observe that if we take the 4-divergence of both sides of 
the inhomogeneous Maxwell equations: d a d^Ff3 a = fi d a J a = the left hand side vanishes 
because again, a symmetric differential operator is contracted with a completely antisym- 
metric field strength tensor. Thus d a J a = 0, which, by some strange coincidence, is the 
charge- current conservation equation in four dimensions. Do you get the feeling that some- 
thing very deep is going on? This is what I love about physics. Beautiful things are really 
beautiful!" I think there is no strange coincidence. This is what I love about physics! The 
beautiful equation d a J a = representing charge conservation is the fundamental symmetry 
behind the inhomogeneous and homogeneous Maxwell equations. "That something very 
deep" is just charge conservation. To consider the continuity equation as the basic axiom 
to build the set of equations known as Maxwell's equations is then a very natural idea. In 
this paper we have followed this idea by obtaining the covariant form of Maxwell's equations 
from the covariant form of the continuity equation. 

It can be argued that other physical assumptions like the use of the Minkowski spacetime 
with its speed of light c and its associated retarded Green function G are ingredients im- 
plicit in the covariant derivation of Maxwell' equations presented here. However, if we first 
postulate the validity of the continuity equation in any four-space, the particular use of the 
Minkowski spacetime (together with its associated retarded Green function G and speed c) 
would not really be a new postulate, but merely an application of our initial postulate. 

We advocate an axiomatic presentation of Maxwell's equations in undergraduate and 
graduate courses of electromagnetism. According to Obradovic [19]: "The axiomatic 
method offers the shortest way to the essence of any theory, enables its more accurate 
formulation and more profound and complete interpretation." In the case of Maxwell's 
equations the essence is charge conservation which should be considered the fundamental 
postulate for an axiomatic presentation of these equations. In the practice, undergraduate 
instructors can find the basis for this presentation in equations (4a) and (4b) of Ref. 2 and 
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graduate instructors in equation (7) of this paper. 
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